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ABSTRACT 

Consider a 2 x 2 factorial experiment with more than 1 replicate. Suppose that we 
have uncertain prior information that the two-factor interaction is zero. We describe 
new simultaneous frequentist confidence intervals for the 4 population cell means, 
with simultaneous confidence coefficient 1 — a, that utilize this prior information 
in the following sense. These simultaneous confidence intervals define a cube with 
expected volume that (a) is relatively small when the two-factor interaction is zero 
and (b) has maximum value that is not too large. Also, these intervals coincide with 
the standard simultaneous confidence intervals obtained by Tukey's method, with 
simultaneous confidence coefficient 1 — a, when the data strongly contradict the prior 
information that the two-factor interaction is zero. We illustrate the application of 
these new simultaneous confidence intervals to a real data set. 
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1. Introduction 

Consider a 2 x 2 factorial experiment with c replicates, where c > 1. Label the 
factors A and B. Suppose that the parameters of interest are the four population cell 
means 9 00 ,9 w ,9 01 ,9u where, for example, 9 W denotes the expected response when 
factor A is high and factor B is low. Also suppose that, on the basis of previous 
experience with similar data sets and/or expert opinion and scientific background, 
we have uncertain prior information that the two-factor interaction is zero. Our aim 
is to find simultaneous frequentist confidence intervals for the population cell means, 
with simultaneous confidence coefficient 1 — a, that utilize this prior information. 

Similarly to Hodges and Lehmann [1], Bickel [2], Farchione and Kabaila [3], 
Kabaila and Giri [4,5] and Kabaila and Tuck [6], our aim is to utilize the uncer- 
tain prior information in the frequentist inference of interest, whilst providing a 
safeguard in case this prior information happens to be incorrect. Throughout this 
paper, we find that the simultaneous confidence intervals of interest define a cube. 
For convenience, we henceforth refer to this cube, rather than the corresponding 
simultaneous confidence intervals. Let 9 = (9 o,9 w ,9 i,9n). The standard 1 — a 
confidence cube for 9 is found using Tukey's method (described e.g. on p. 289 of [7]). 
We assess a 1 — a confidence cube for 9 using the ratio (expected volume of this 
confidence cube) /(expected volume of standard 1 — a confidence cube). We call this 
ratio the scaled expected volume of this confidence cube. We say that this confi- 
dence cube utilizes the prior information if it has the following desirable properties. 
This confidence cube has scaled expected volume that (a) is significantly less than 1 
when the two-factor interaction is zero and (b) has a maximum value that is not too 
much larger than 1. Also, this confidence cube coincides with the standard 1 — a 
confidence cube when the data strongly contradict the prior information that the 
two-factor interaction is zero. 

An attempt to utilize the uncertain prior information is as follows. We carry out a 
preliminary test of the null hypothesis that the two-factor interaction is zero against 
the alternative hypothesis that it is non-zero. If this null hypothesis is accepted then 
the confidence cube for 9, with nominal confidence coefficient 1 — a, is constructed 
assuming that it is known a priori that the two-factor interaction is zero; otherwise 
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the standard 1 — a confidence cube is used. We call this the naive 1 — a confidence 
cube for 9. This assumption is false and it leads to a naive 1 — a confidence cube 
with minimum coverage probability less than 1 — a. For example, for a = 0.05, 
c = 2 and a preliminary test with level of significance 0.05, this minimum coverage 
probability is 0.9078. The poor coverage properties of these naive confidence cubes 
are presaged by the following two strands of literature. The first strand concerns the 
poor properties of inferences about main effects after preliminary hypothesis tests 
in factorial experiments, see e.g. Neyman [8], Traxler [9], Bohrer and Sheft [10], 
Fabian [11], Shaffer [12] and Ng [13]. The second strand concerns the poor coverage 
properties of naive (non-simultaneous) confidence intervals in the context of linear 
regression models with zero-mean normal errors, see e.g. Kabaila [14, 15], Kabaila 
and Leeb [16], Giri and Kabaila [17] and Kabaila and Giri [18]. 

Whilst the naive 1 — a confidence cube fails to properly utilize the prior informa- 
tion, its form (described in Section 2) will be used to provide some motivation for 
the new 1 — a confidence cube that utilizes the uncertain prior information and is 
described in Section 3. In this section we also provide a numerical illustration of the 
properties of this new confidence cube for 1 — a = 0.95 and c = 2. The two-factor 
interaction is described by the parameter in the regression model used for the 
experiment. The uncertain prior information is that /3 12 = 0. Define the parameter 
7 = P12 /\j var(/3 12 ), where /3 12 denotes the least squares estimator of (3 l2 . As proved 
in Section 3, the scaled expected volume of the new confidence cube for 9 is an even 
function of 7. The bottom panel of Figure 2 is a plot of the square root of the 
scaled expected volume of the new 0.95 confidence cube for 9, as a function of 7. 
When the prior information is correct (i.e. 7 = 0), we gain since the square root 
of the scaled expected volume is substantially smaller than 1. The maximum value 
of the square root of the scaled expected volume is not too large. The new 0.95 
confidence interval for 9 coincides with the standard 1 — a confidence cube when the 
data strongly contradicts the prior information. This is reflected in Figure 2 by the 
fact that the square root of the scaled expected volume approaches 1 as 7 — * 00. In 
Section 4 we illustrate the application of the new 1 — a confidence cube to a real 
data set. 
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2. The naive 1 — a confidence cube 

Let Y denote the response and x 1 and x 2 denote the coded levels for factor A 
and factor B respectively, where x\ takes values —1 and 1 when the factor A takes 
the values low and high respectively and x 2 takes values —1 and 1 when the factor 
B takes the values low and high respectively. We assume the model 

Y = p + Pix ± + p 2 x 2 + p 12 x ± x 2 + e (1) 

where (3 , Pi, /3 2 and Pi 2 are unknown parameters and the e for different response 
measurements are independent and identically N(0, a 2 ) distributed. Because we are 
considering c replicates, the number of measurements of the response is n = 4c. 
The dimension of the regression parameter vector (/3 , Pi, p 2 , Pu) is p = 4. The 
parameters of interest are 9 00 = Po — Pi — P2 + Pi 2 , W = Po + Pi — P 2 — Pi 2 , 
$01 — Po — Pi + P2 — P12 and On = Po + Pi + P2 + Pu- The uncertain prior information 
is that Pi 2 = 0. 

The naive 1 — a confidence cube is defined as follows. We carry out a preliminary 
test of the null hypothesis that Pi 2 = against the alternative hypothesis that 
P12 7^ 0. If this null hypothesis is accepted then the confidence cube is constructed 
assuming that it was known a priori that Pi 2 = 0; otherwise the standard 1 — a 
confidence cube for is used. 

Let 0o, Pi, P 2 ,Pi 2 ) denote the least squares estimator of {Po, Pi, P2, Pu)- Also, 
let (6 o, 6ioj @oi) @n) denote the least squares estimator of (0 00 , W , O i, On) . Note 
that (000) @ioj @oij ©11,^12) has a multivariate normal distribution with mean 
(0oo,0io,0oi,0n, Pi 2 ) and covariance matrix a 2 V, where 

'4 1 " 

4 0-1 
4 0-1 
4 1 
1-1-11 1 

Let Vij denote the (i,j) th element of V. Also, let a 1 denote the usual estimator of 
a 2 obtained by fitting the the full model to the data. Define W = a jo. Note that 
W has the same distribution as \jQj{n — p) where Q ~ Xn-p- Let fw denote the 
probability density function of W. Define di- a by P(maxj = i ; ... )4 \Zi\/W < di- a ) = 

4 



V 
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1 — a, where Z\, Z2, Z 3 , Z4 and W independent random variables and Z^ ~ N(0, 1) 
(i = 1,2,3,4). We use [a ± b] to denote the interval [a — b,a + b]. The standard 
1 — a confidence cube for 9 is 



0oo ± di-aJv^a 



x 



©io ± dx-ay/v^a 

601 ± d^Jv^a 



x 



X 



0ii ± d\- a Jv^a 



The naive 1 — a confidence cube for 9 is obtained as follows. The usual test 
statistic for testing the null hypothesis that f3\2 = against the alternative hy- 
pothesis that P12 7^ is P12/ (d'\^V55)- This test statistic has a t„_ p distribution 
under this null hypothesis. Suppose that, for some given positive number q, we 
fix P12 at if |/3i2|/(<5"v / ^55) — ?! otherwise we allow P± 2 to vary freely. Define 
a 2 = ((n — p)cr 2 + 2 2 /v 55 )) /(n — p + 1). This is the usual estimator of a 2 ob- 
tained by fitting the full model to the data when it is assumed that [3\2 = 0. Define 
W = a /a. Note that W has the same distribution as \J Q/(n — p + 1), where 
Q ~ Xn-p+i when p u = 0. Suppose that Z ± , Z 2 , Z 3 and W are independent 
random variables and Z { ~ iV(0, 1) (i = 1,2,3). Let Z 1 = (Z 1 - Z 2 - Z 3 )/y/3, 
Z 2 = {Z 1 + Z 2 - Z 3 )/y/3, Z 3 = (Z 1 -Z 2 + Z 3 )/V3 and Z 4 = (Z 1 + Z 2 + Z 3 )/V3. 
Define d\_ a by P(max i= i j .. .,4 < di_ a ) = 1 — a. The naive 1 — a confidence 

cube is obtained as follows. If 0\2\l{p\fv^) > Q then this confidence cube is /. If, 
on the other hand, \f3\2\/{ffyjvsi) < Q then this confidence cube is 

[0oo - P12 ± ( VsTc/2)di-ad-] x [0io + /3i 2 ± (^Tc/^d^a] 

x [0oi + $ 12 ± (v / 3A/2)di_ a a] x [0no - /3i2 ± (v / 37c/2)rfi-«^] 
The naive 1 — a confidence cube can be expressed in the form 

P12 \ , ^ IA2I 







00 



x 



1>11<7 6 



± \/Vu& s 



©10 + y^crb 



P 



12 



X 



©01 + v^ 5 " & 



012 
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± V^44 ° S 







12 



(2) 
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where 




for \x\ > q 

x/2 for \x\ < q. 

di- a for x > q 



di- a (y/3/2)y/(n -p + x 2 )/(n-p + 1) for < x < q. 



3. New 1 — a confidence cube utilizing the prior information 

We introduce a confidence cube for 9 that is similar in form to the naive 1 — a 
confidence cube, described in the previous section, but with a great "loosening up" 
of the forms that the functions b and s can take. Define the confidence cube J(b, s) 
for 9 to have the form (j2J), where the functions b and s are required to satisfy the 
following restriction. 

Restriction 1 b : R — > K is an odd function and s : [0, oo) — > [0, oo). 

Invariance arguments, of the type used by Farchione and Kabaila [3] and Kabaila 

and Giri [4] , may be used to motivate this restriction. For the sake of brevity, these 

arguments are omitted. We also require the functions b and s to satisfy the following 

restriction. 

Restriction 2 b and s are continuous functions. 

This implies that the endpoints of the simultaneous confidence intervals correspond- 
ing to the confidence cube J(b, s) are continuous functions of the data. Finally, 
we require the confidence cube J(b, s) to coincide with the standard 1 — a confi- 
dence cube I when the data strongly contradict the prior information. The statistic 
|^i2|/(<3"a/^22) provides some indication of how far away Pu/ (cr^/v^) is from 0. We 
therefore require that the functions b and s satisfy the following restriction. 

Restriction 3 b(x) = for all |x| > r and s(x) = d\_ a for all x > r where r is a 
(sufficiently large) specified positive number. 



Define 7 = /W^v 7 ^ T i = (A) ~ A))/(<V^ T 2 = 0i - Pi)/(<Ty/v£), 
T ?J = ((3 2 - /5 2 )/(o- v /: ^5) and H = /WOv 7 ^)- Tnen define G ± = {T x - T 2 - 
T 3 + H- 7 )/2, G 2 = (T 1 + T 2 -T 3 -H + 7 )/2, G 3 = (T a - T 2 + T 3 - H + 7 )/2, 
G 4 = (Tx + T 2 + T 3 + H - 7)/2. Note that H ~ A^( 7 , 1). It is straightforward to 
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show that the coverage probability P{6 e J(b, s)) is equal to 



P(£(H, W) < d < W), W) < -G 2 < u(H, W), 

£(H, W) < -G 3 < u(H, W), £(H, W) < G A < u(H, W)) (3) 

where the functions £(■, •) : R x [0, oo) — > R and «(•, •) : R x [0, oo) — > R are defined 
by £(h, w) = b{h / 'w)w — s(\h\ / 'w)w and u(h, w) = bihj w)w + s(\h\ / w)w . As we prove 
later, for given b and s, the coverage probability of J{b, s) is an even function of 7. 
We denote this coverage probability by 0(7; b, s). 

Part of our evaluation of the confidence cube J(b, s) consists of comparing it 
with the standard 1 — a confidence cube / using the criterion 

expected volume of J{b, s) 
expected volume of / 

We call this the scaled expected volume of J(b, s). This is equal to 

di- a E(W*) ■ 

This is a function of 7 for given s. We denote this function by e(7; s). Clearly, for 
given function s, e(7; s) is an even function of 7. When J(6, s) is a 1 — a confidence 
cube, we use \/e( / y; s) to measure the efficiency of the standard 1 — a confidence 
cube relative to J(b, s), for parameter value 7. The reason for this is that a/ e{j; s) 
is a measure of the ratio of sample sizes required for the expected volumes of these 
two confidence cubes to be equal. 

Our aim is to find functions b and s that satisfy Restrictions 1-3 and such that 
(a) the minimum of 0(7; 6, s) over 7 is 1 — a and (b) the weighted average 

/oo 
(e( 7 ; s) - 1) <Mt) (4) 
-00 

is minimized, where the weight function v has been chosen to be 

v(x) = Xx + H(x) for all i£l, (5) 

where A is a specified nonnegative number and 7i is the unit step function defined by 
H(x) = for x < and H(x) = 1 for x > 0. The larger the value of A, the smaller 
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the relative weight given to minimizing e(7; s) for 7 = 0, as opposed to minimizing 
e(7; s) for other values of 7. The idea of minimizing a weighted average expected 
length of a confidence interval, subject to a coverage probability constraint, is due 
to Pratt [19]. The particular weight function (JHJ) was first used in related contexts 
by Farchione and Kabaila [3] and Kabaila and Giri [4]. 

The following theorem provides computationally convenient expressions for the 
coverage probability and scaled expected length of Jib, s). 

Theorem 1. 

(a) Define li = 2£(h,w)+t 3 -h+>y, u x = 2u(h,w)+t 3 -h+>y, £ 2 = -2u(h,w)+t 3 +h- 
7 and u 2 = —2£(h, w) + t 3 + h — 7. Now define t\ = max(^i, — U\), u\ = min(iti, — ft), 
£2 = max(£ 2 , — u 2 ) and u 2 = min(u 2 , —t<i). We use these functions to define 



k(t 3 ,h,w,i) 



if either £1 > U\ or £ 2 > u 2 , 
$(u x /\/2) -$(4/v^)) ($(u 2 /V2) -$(£ 2 /V2)^ otherwise, 



where $ denotes the N(0, 1) distribution function. 

Define ft (w) = —di^ a w andv){w) = d\- a w. Also define £\ = 2ft (w) +t 3 — ft. + 7, 
u\ = 2v){w) +t 3 -ft + 7, 4 = -2w t (w) + t 3 + ft-7 and u\ = -2ft (w) + 1 3 + ft - 7. 
Now define I\ = max(4, — u\), u{ = mm(u\, —£{) , £ 2 = max(4, — u 2 ) and u 2 = 
min(w2, -4)- We use these functions to de&ne 

if either l\ > u\ or l\> u\, 

^(u\/V2) - $(4/V2)) (<S>(u\/\/2) - $(4/V2)) otherwise. 



k ] {t 3 ,h,w,^) 



The coverage probability of J{b, s), denoted by 0(7; b, s), is equal to 



00 



(l—a)+ / / (k(t 3 ,wx,w,^)—k\t 3 ,wx,w,^))(f)(t 3 )dt 3 (j)(wx—'y)dxwfw('w)dw 

J J —v J —00 

(6) 

where (f) denotes the N(0, 1) probability density function. For given functions b and 
s, 0(7; b, s) is an even function of 7. 

(b) The scaled expected volume of J(b, s), denoted by e(7; s), is equal to 

poo rr 

l + J\ — E(W 4 ) J J <t>{wx-l)dxw 5 f w (w)dw. (7) 
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This theorem is proved in the Appendix. Substituting §7§ into (J4]), we obtain that 
(jl]) is equal to 

1 



1 



oo JO 



(s (\x\) — d*_ a ) 4>{wx — 7) dxw 5 fw( w ) dwdu^j) 



df_ a E(W*) 
2 



(AH) -<£-«) / 

J — c 



(■ira — 7) £^(7) dxw 5 fw( w ) dw 



(s (x) — d^_ a ) (A + cf)(wx)) dxw 5 fw{w) dw 



For computational feasibility, we specify the following parametric forms for the 
functions 6 and s. We require 6 to be a continuous function and so it is necessary that 
6(0) = 0. Suppose that x±, . . . ,x m satisfy = x\ < 22 < • • • < x m = r. Obviously, 
b(xi) = 0, b(x m ) = and s(x m ) = di- a . The function b is fully specified by the 
vector (6(0:2), • • • > b(x m -i)) as follows. Because 6 is assumed to be an odd function, 
we know that Xj) = — 6(xj) for % = 2, . . . ,m. We specify the value of b(x) for 
any x G [— r, r] by cubic spline interpolation for these given function values, subject 
to the constraint that 6'(— r) = and 6'(r) = 0. We fully specify the function s by 
the vector (s(xi), . . . , s(x m _i)J as follows. The value of s(x) for any x G [0, r] is 
specified by cubic spline interpolation for these given function values (without any 
endpoint conditions on the first derivative of s). We call x±,X2, ■ ■ -x m the knots. 

To conclude, the new 1 — a confidence cube for 9 that utilizes the uncertain 
prior information that /3 12 = is obtained as follows. For a judiciously-chosen set 
of values of r, A and knots x iy we carry out the following computational procedure. 

Computational Procedure 

Compute the functions 6 and s, satisfying Restrictions 1-3 and taking the parametric 
forms described above, such that (a) the minimum over 7 > of (jSJ) is 1 — a and 



(b) the criterion (jSj) is minimized. Plot y~e(pi\ s) as a function of 7 > 0. 

Based on these plots and the strength of our prior information that /? 12 = 0, we 
choose appropriate values of r, A and knots x,- t . The confidence cube corresponding 
to this choice is the new 1 — a confidence cube for 9. 

Consider the case that c = 2 and 1 — a = 0.95. We followed this Computational 
Procedure, with r = 8, A = 0.08 and evenly-spaced knots Xi at 0,r/6, . . . ,r. The 
resulting functions 6 and s, which specify the new 0.95 confidence cube for 9 are 
plotted in Figure 1. The performance of this confidence cube is shown in Figure 
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2. This confidence cube has the attractive property that its coverage probability is 
0.95 throughout the parameter space. When the prior information is correct (i.e. 



7 = 0), we gain since ^/e(0; s) = 0.8558. The maximum value of ^/ e (7! s ) is 1-0956. 
This confidence cube coincides with the standard 0.95 confidence cube for when 
the data strongly contradicts the prior information, so that a/ e(7; s) approaches 1 
as 7 — > oo. 



b(x) 




s(x) 




Figure 1: Plots of the functions b and s that specify the new 0.95 confidence cube 
for 9 = (9 00 , 0io, 0oi, 012), when c = 2, r = 8 and A = 0.08. 
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Figure 2: Plots of the coverage probability and the square root of the scaled expected 

volume e(7; s) (as functions of 7 = (3u/ y var(/3i2)) of the new 0.95 confidence cube 
for 6 = (9 0Q , 810, 9 i, 6i 2 ). This cube was obtained using c = 2, r = 8 and A = 0.08. 
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4. Illustration of the application of the new confidence cube 

In this section we illustrate the application of the new 1 — a confidence cube, 
utilizing the uncertain prior information, to a real data set. We extract a 2 x 2 
factorial data set from the 2 3 factorial data set described in Table 7.5 of Box et al 
[20] as follows. Define x\ — — 1 and X\ = 1 for "Time of addition of HNO3" equal to 
2 hours and 7 hours, respectively. Also define x 2 = — 1 and X2 = 1 for "heel absent" 
and "heel present", respectively. The observed responses are the following: 



We use the model ([T]). The discussion on p. 265 of Box et al [20] implies that there is 
uncertain prior information that ft\ 2 = 0. The discussion on p. 266 of Box et al [20] 
implies that there is an estimator a 2 of a 2 , obtained from other related experiments, 
with the property that a 2 jo 2 ~ Q/m where Q ~ xin an d m is effectively infinite. 
The observed value of a is 0.8. The standard 0.95 confidence cube for 6 is 



[87.2 ± 1.99272] x [88.4 ± 1.99272] x [86.7 ± 1.99272] x [89.2 ± 1.99272]. 

We have also computed the new 0.95 confidence cube for 6, using d — 6, A = 0.08 
and equally-spaced knots at 0, 1, . . . , 6. This new confidence cube is 

[87.1748 ±1.88504] x [88.4252 ± 1.88504] x [86.7252 ± 1.88504] x [89.1748 ± 1.88504]. 

Thus 



For this data set, we have clearly gained by using the new 0.95 confidence cube for 
0. 

Remark The new 1 — a confidence cube is computed to satisfy the constraint that 
its minimum coverage probability is 1 — a. For the examples described in both the 
previous section (for which c = 2, r = 8 and A = 0.08) and the current section 
(for which c is effectively infinite, r = 6 and A = 0.08), it is remarkable that the 



For (xi, x 2 ) 
For (xi ) x 2 ) 
For (xx : x 2 ) 
For (x 1 ,x 2 ) 



(-1,-1), y = 87.2 
(1,-1), y = 88.4. 
(-1,1), 3/ = 86.7. 
(1,1), y = 89.2. 
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new 0.95 confidence cube has coverage probability equal to 0.95 throughout the 
parameter space. The new 0.95 confidence cube was also computed for (a) c = 3, 
r = 8 and A = 0.08, (b) c = 7, r = 8 and A = 0.08 and (c) c = 20, r = 8 and 
A = 0.08. In each of these cases, the new 0.95 confidence cube also has coverage 
probability equal to 0.95 throughout the parameter space. This provides strong 
empirical evidence that the new 1 — a confidence cube has the attractive property 
that its coverage probability is equal to 1 — a throughout the parameter space. 

Appendix. Proof of Theorem 1 

In this appendix we prove Theorem 1. 
Proof of part (a). 

Remember that 0(7; b, s) is equal to ([3]). It follows from the N(j, 1) distribution of 
H and the independence of the random vectors (G x , G 2 , G 3 , G4, H) and W that (J3j) 
is equal to 



a(h,w) =P(£(h,w) < G x < u(h,w), £(h,w) < -G 2 < u(h,w), 

£(h, w) < -G 3 < u{h } w), £(h, w)<G A < u(h, w)\H = ti). 

Recall the expressions for G\, G2, G3 and G4 in terms of T\, T2, T3 and H. Also 
note that 7\, T 2 , T 3 and H are independent random variables and that I\, T 2 and 
T 3 are identically A^(0, 1) distributed. Thus T\ — T 2 and T\ + T 2 are independent 
A^O, 2) distributed random variables and a(h, w) is equal to 



i P{t x < 21 - T 2 < u h -«i < Ti - T 2 < -£ 1} 

—00 

£2 < + T 2 < u 2 , -u 2 <T, + T 2 < -l 2 ) <f>(t 3 ) dt 3 




where 





00 



P{ii < T x - T 2 < u x )P(t 2 < T x + T 2 < u 2 ) 0(t 3 ) dt 3 



—00 



00 




—00 



Thus 0(7; b, s) is equal to 




(A.l) 
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The standard 1 — a confidence cube I has coverage probability 1 — a. Hence 



OO /"OO /"OO 



l-a= / / k j {t 3 ,h,w,^)(p{t3)dt 3 (j){h-^)dhf w {w)dw. (A.2) 

Jo J — oo J — oo 

Subtracting ( 1A.2I) from ( 1A.1I) and noting that, by Restriction 3, b(x) = for all 
\x\ > r and s(x) = for all x > r, we find that 0(7; b, s) is equal to 

/»oo /*rtL> /»oo 

(!-«)+/ / / (k(t 3 ,h,w^)-k\t 3 ,h,w-f)) <j)(t 3 )dt 3 (f)(h-^)dhf w (w)dw. 



' •/ — rui »/ —00 

Changing the variable of integration from h to x = h/w, we obtain ([6]). It is 
straightforward to show that P{9 e J(6, s)) is an even function of 7. 

Proof of part (b). 

The random variables H and W are independent. Thus 

e( 7 ; a) = ^ jj^FIj /" s * ) W " 7) * ^ 4 M«0 ^ (A.3) 
Obviously, 

OO /* OO 

4 jl/l .A JL...4 







1 = df E{w4 , I I d\_ a <j ) {h- 1 )dhw^ w {w)dw. (A.4) 



Note that s(x) = d±^ a for all x > r. Subtracting (I A .40 from (1A.3j) . we therefore 
obtain 

e(7;s) = '^fjC ( s4 («) " ^-^ 4 ^m 

Changing the variable of integration from ft, to x = h/w, we obtain (171) . 
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